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We study the conditions for the existence of unpaired Majorana modes at the ends of vortex lines or 
the side edges of a layered topological superconductor. We show that the problem is mapped to that 
of a general Majorana chain and extend Kitaev's condition for the existence of its nontrivial phase by 
providing an additional condition when a supercurrent flows in the chain. Unpaired Majorana bound 
states may exist in a vortex line that threads the layers if the spin-orbit coupling has certain in-layer 
components but, interestingly, only if a nonzero supercurrent is maintained along the vortex. We 
discuss the exchange statistics of vortices in the presence of unpaired Majorana modes and comment 
on their experimental detection. 



I. INTRODUCTION 

Majorana fermions were introduced in 1937^ as neu- 
tral spin-1/2 particles that are their own anti-particles. 
Mathematically, they are described by a modified two- 
component Dirac equation, which admits purely real 
solutions. Majorana himself speculated that neutrinos 
might be described by this formalism. Indeed, after many 
years of hiatus, this appears now to be a promising sug- 
gestiorP. Other candidate Majorana fermions are found 
in supersymmetric theories and theories of dark mat- 
tei^^l. However, the experimental observation of Majorana 
fermions in these theories seems distant. 

In recent years, an increasing number of proposals have 
been made for the realization of Majorana fermions as 
quasiparticles in condensed matter system^^HIl, While 
these proposals vary in their implementation, they gener- 
ically exploit the special structure of quasiparticles in a 
superconductor, each of which is a mixture of an electron 
and a hole. By virtue of this structure, a quasiparticle 
whose energy is in the middle of the superconducting gap 
(a "zero mode" pinned to the chemical potential) is au- 
tomatically a Majorana fermion. Hence, the aim of these 
proposals is to find systems where such zero modes can 
be stabilized. The existence of zero modes is a topologi- 
cal property of the Hamiltonian: they are robust against 
smooth changes in potentials. This topological protec- 
tion would be an advantage in the experimental search for 
Majorana fermions once the right Hamiltonian is found 
or engineered. 

A prime example is provided by a chiral p-wave {px ± 
ipy) superconductor in two dimension^. A thin layer of 
the A-phase of superfluid ^He has this pairing symmetry. 
There is also evidence that Sr2Ru04, a layered material, 
may realize a similar pairing state in each layei^. A close 
relative of this state, which combines both chiralities into 
a time-reversal invariant pairing state, also hosts Ma- 
jorana modes protected by the additional time-reversal 
symmetry. The B-phase of superfluid "^He realizes this 
pairing symmetry. The doping-induced sup erco nducting 
state of the topological insulator Cua;Bi2Se5}2li2Janother 
layered material, might be another realization^. Vor- 
tices in these topological superconductors support Majo- 
rana bound state^^. In addition, the boundaries of these 




FIG. I. (color online) The experimental setup with a mul- 
tilayered topological superconducting annulus (green). The 
hole is threaded by the magnetic field B shown by the (or- 
ange) straight arrow. Superconducting leads, shown by (yel- 
low) bars are attached to the top and bottom surface. A su- 
percurrent js shown by (blue) curvy arrow, flowing along the 
hole, can be maintained between the leads. Majorana edge 
states are marked by arrows to show the direction of propa- 
gation. In the topological state, unpaired Majorana fermions 
are bound to edges at the top and bottom layer, shown by 
thick black circles. 



systems also contain propagating Majorana modes. 

In two spatial dimensions, the ground state of the 
superconductor with several Majorana bound states is 
highly degenerate and is unitarily mapped within the 
degenerate subspace upon exchange of two Majorana 
bound states. Therefore, these Majorana bound states 
behave as non-Abelian anyoni^^ and could be useful for 
fault-tolerant quantum computatior^. However, these 
properties may not survive in a layered system. Indeed, 
while two spatially separated Majorana fermions are pro- 
tected, once they reach tunneling proximity they may 
hybridize into a regular fermion. So, if tunneling ampli- 
tude between Majorana fermions in and between layers is 
nonzero it might be expected that no unpaired Majorana 
fermions would survive. 

In this paper, we study the conditions under which 
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unpaired Majorana fermions exist in layered topological 
superconductors. Specifically, we study the Majorana 
fermions in a vortex line or the boundaries of the system 
normal to the layers (see Fig. [T]) creating a "Majorana 
chain" along the layers. In a spin-triplet chiral p-wave su- 
perconductor, a stable half-quantum vortex (HQV) car- 
rying a flux h/Ae is allowed which supports a single Ma- 
jorana bound state in each layer. Recently, strong ev- 
idence for the existence of HQVs in Sr2Ru04 has been 
reported^. A full-quantum vortex (FQV) carrying a flux 
/i/2e, on the other hand, supports two Majorana bound 
states in each layer. The vortex in a time-reversal in- 
variant topological superconductor also carries h/2e flux 
and two Majorana bound states. The boundaries of both 
systems host two propagating Majorana modes per layer. 
Motivated by these observations and by a recent proposal 
to exploit h/2e FQVi^, we focus on the problem of two 
Majorana fermions per layer. For brevity we limit our- 
selves to the case of a chiral p-wave superconductor but 
our study may be readily generalized to other situations. 

Our main finding is that passing a supercurrent in the 
direction normal to the layers could stabilize unpaired 
Majorana bound states at the surface layers of the sys- 
tem. The equivalent Majorana chain contains complex 
fermion tunneling amplitudes. Therefore, we generalize 
Kitaev's original work'' to this case. Our result can be 
stated as follows. The intra-layer Majorana tunneling 
tends to pair the two Majorana fermions in each layer 
into a regular fermion and thus acts as a chemical po- 
tential. The interlayer tunneling amplitude, on the other 
hand, tends to pair Majorana fermions in different layers 
and acts as both complex tunneling as well as Cooper 
pairing of the regular fermions. The Majorana chain is 
then equivalent to a fermion chain in a superconduct- 
ing state with a steady supercurrent. The condition for 
unpaired Majorana fermions at the chain boundaries is 
that the chemical potential lies in the fermion bandwidth 
(same as the condition derived by Kitaev) and that the 
supercurrent is below the critical current of the supercon- 
ducting state of the fermion chain (in addition to Kitaev's 
condition) . 

In a layered chiral p-wave superconductor we find that, 
remarkably, unpaired Majorana fermions may only exist 
when a supercurrent is maintained normal to the layers. 
We shall discuss a set of sufficient conditions to stabi- 
lize such unpaired Majorana fermions. We explore their 
non-Abelian exchange statistics and also discuss ways to 
detect them. Details of our calculations are presented in 
two appendices. 



II. THE MAJORANA CHAIN 

Let us start with the following problem: consider a 
one-dimensional lattice with sites t — 1,2^ ■ ■ ■ N and two 
Majorana fermions 7crf , cr S {f, i} at each site. The most 
general Hamiltonian describing Majorana tunneling at 




FIG. 2. The sign assignments 7o-f — ^ ±7(t<, f" £ {tii} of M-a.- 
jorana operators for Z2 spectral symmetries. Black and white 
circles indicate opposite sign assignments. Dashed (solid) 
lines show tunneling amplitudes that are multiplied with a — 
(+) sign, (a) — >■ —iJ,e,te t}, A*; — >■ A|; and (b) tt — > —te, 



each site and between nearest neighbors is 

N-l N 
e=l a-a-' 1=1 

with tca'iTCti G The Majorana operators satisfy the 
algebra {'ycri,l<j'i'} — 2Sa-a-'Su'- Let us define a fermionic 
operator for each layer, ipg — ^{'jfe + ijii) which sat- 
isfy the usual fermionic algebra {ipl , ipg' } = 6u' and 
{ipi,ipt} = 0. Then 



1=1 1=1 

(2) 

where 

ti = t^ii ~ h^i + iittn + hu) = \tiW''' , (3a) 
Ai = t^ii + t^^i + im^i - t^ii) = \Ai\e"^\ (3b) 
m = -2ae, (3c) 

and c = — J2eLi is a constant. 

The spectral properties of Hamiltonians ([l]) and ([2]) 
are invariant under the following operations: (a) /i£ — > 
-im^ti t\, At ^ A\; and (b) tg -ti, Ai -Ai. 
In each case, the changes can be undone by a redefinition 
of the Majorana operators •j^re — ±7(t£- The choice of 
signs is illustrated in Fig. [2j 

In general, we can remove the complex phase of ti by 
a gauge transformation 



(4a) 



A, -> A,e^(''^+''^+i) .-. 0, ^ 0, -h {m + %+i)- (4b) 

It follows from the last equation that there is a 0^ -|- 
^^_i_i(mod 27r) contribution to the interlayer supercurrent, 
js/ oc (pi+i ~ 4'ii which is nonzero when ti ^ M. 

From here on let us assume that all amplitudes are 
constant, = t,Ai — A,/i£ — fi. Then, the spectrum 
of ([2]) as a function of the lattice momentum q is 



Eq± ^ (Imt)sinq ± ^ [(Ret) cos g - /i/2]"' + jApsin^ q, 

(5) 
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which has the symmetry E- 



-E, 



We can see 



from this expression that the Im (t) shifts the energies 
and therefore gives rise to a current. 

The Hamiltonian in (2]) when t € M was considered by 
Kitaev^, who showed that it admits two phases. In the 
trivial phase, there are no unpaired Majorana modes. In 
the nontrivial phase, there are unpaired Majorana modes 
at the boundaries of the system. The nontrivial phase 
exists wherP 



\n\ < 2\Ret\ 



(6) 



which is the condition for a partially filled band. 

We now show that when t ^ M., the nontrivial phase 
obtains when in addition to (|6|, 



lAI > llmtl. 



(7) 



The quantity |A| — \lmt\ is proportional to the gap at 
the Fermi points where (Ret) cos qp — fi/2. Physically, 
this condition ensures that the supercurrent is less than 
the critical value above which superconducting gap in ^ 
is destroyed: the complex phase of t is restricted to 
I tan 1 < 2|A|/|/^|, which guarantees that the supercur- 
rent js oc 29 (for small | A|) is less than the critical current 
in the chairPSJ. This is one of our main results. 

To see this, we pass to the continuum limit of ([T]), 
by writing z = ia, a — 0. There are two counter- 
propagating low-energy modes with lattice momenta 
g sa and q k, ja where the gap vanishes. Expand- 
ing around their respective momenta we have 7(t£+i ~ 
±[7,(z)+5,7,(z)]. So, 

H,= '-Jdz [^^ {Td, + M) Co + C (-ra. + M^) ^ , 



r = {lnit)I + (Im A)(T^ + (Re A)cr^ 



M 



Ret - -fi] lay, 



(8a) 
(8b) 

-Ret — -/i j i(Tj,,(8c) 



with ^o(7r) = (7t;7^)"^ fo'" modes near (tt). We now 
look for zero-energy solutions of the Dirac operator D = 
Tdz + M localized where the mass m = Re t — changes 
sign. The localized mode corresponds to an unpaired 
Majorana bound state at the boundary of the nontrivial 
phase. Let us take m = |TO|sgn(z) for simplicity; then, if 



eEE |Ap-(Imt)' >0, 
we find ZJe"'^'^'^ = with k = |m|/^/e, and 



(9) 



sgn(m)v/e- 
Im(A- 



- Re A 

t) 



Im (A - t) 
sgn(m)Y^ -I- Re A 



(10) 

whichever is nonzero. For modes near tt the Dirac opera- 
tor is = —Ydz + M^. So, for the mass changes 
sign either for the modes near or the ones near tt/o and 
there is only one zero mode. Since the sign of t and /i 
can be changed by appropriate sign assignments of the 
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FIG. 3. (color online) The spectrum and eigenstates for a 
Majorana chain with N — 100 sites. Top: The spectrum of 
the Majorana chain as function of x = tfi (not all eigenvalues 
are shown). The other parameters are fixed cLS (y. — 0.2, — 
til = 0.4, tif = 0.7. The conditions Q and ^ yield x < 
0.5 and x > 0.1, respectively. Two zero modes are clearly 
visible in the spectrum for 0.1 < a; < 0.5. Bottom: The wave 
functions of the two zero modes for x — 0.25. Squares (blue) 
and circles (red) distinguish the two modes. The Majorana 
modes are localized at the two ends of the chain. 



Majorana operators, we find that ^ gives the condition 
for the nontrivial phase also when t ^ M. subject to the 
additional condition ([t]). 

We have also confirmed these results numerically. In 
Fig. [Sj we present the results of exact diagonalization of 
Hamiltonian ([I]) on a chain with N — 100 sites. The pa- 
rameters are chosen to explore the conditions ([6| and ([t]) 
for the presence of Majorana modes. It is clear that both 
conditions are necessary (and sufficient) for the existence 
of unpaired Majorana modes. 

We now turn to study in detail the case of a chiral 
p-wave superconductor of the type that is thought to 
be realized in Sr2Ru04. In this state, while the time- 
reversal symmetry is spontaneously broken by the order 
parameter, the presence of spin degeneracy gives rise to 
two Majorana modes per layer. These Majorana modes 
may then couple both within a layer and in between lay- 
ers through the Zeeman term, the spin-orbit term, and 
direct wave function overlaps between layers. We show 



that in a finite range of parameters a topological state is 
established in the vortex line (see Fig. 111). 



III. UNPAIRED MAJORANA BOUND STATES 
IN A VORTEX LINE 

The quasiparticles of the superconducting state are de- 
scribed by the Bogolubov-de Gennes (BdG) Hamiltonian 
^^i?\['/2, where the Nambu operator 'I' = {ip , 'ip'^ , ip is 
the electron annihilation operator, 



iof 



H 



h A 

At -/i^ 



(11) 



with h the Hamiltonian density of the system in the 
absence of superconductivity, and A the pairing order 
parameter of the condensate. We have suppressed the 
spin index for brevity. The order parameter is skew- 
symmetric A^ — —A, therefore TxH*Tx = —H where 
is the Pauli matrix acting in the Nambu space. The BdG 
quasiparticle of energy E is annihilated by = ^ Xe^ ^ 
where Hxe = Exe- By the symmetry of the BdG 
Hamiltonian, X-e — TxXe- Noting that by construc- 
tion = '^'^Tx, we have = ^'-_e. In particular, at 
zero energy, = ^'J = 7 would be a Majorana fermion. 

For a spin-triplet chiral p-wave superconductor the 
pairing term has the form 

A = ^«Ae'*/'{a, - idy, {a ■ d)ay}e'^'\ (12) 

where $ is the phase of the order parameter, and ua is 
a constant gap velocity. Here d determines the direction 
of pairing in the space spanned by the spin-triplet states, 
and {•,•} is the anticommutator. Restoring the spin 
indices, now ^p = {ip^,ip\)'^ . 

The Hamiltonian is symmetric under the mapping 
d — d and <f> $ + tt. This Z2 symmetry stabilizes a 
half-quantum vortex. When encircling a HQV, the phase 
of the order parameter changes by tt, so $ = -I- 6*72 
where (p is Hic smooth part of the phase. The extra an- 
gular momentum needed to produce a single- valued wave 
function is provided by rotating d to — d. For simplic- 
ity, we consider a configuration where d rotates in the 
xy-plane, say d = (cos 0/2, — sin 0/2, 0). For this con- 
figuration, there is a 27r winding of the oder parameter 
phase only in the spin-f block of the BdG Hamiltonian. 
Therefore, there is a single zero mode with support in 
spin-f component in each layer and an unpaired Majo- 
rana bound state, 7^ (see Appendix [A| for details). 

In the full quantum vortex d is constant and the phase 
of the order parameter winds by 2-k. Both spin compo- 
nents now see the phase winding and there are, therefore, 
two zero modes per layer corresponding to two Majorana 
bound states, 7^ and 7^. 

The Majorana fermions can tunnel between the lay- 
ers and, in the case of the FQV, within each layer. 
Here we consider first-order perturbations that break 
spin rotation symmetry, namely the Zeeman splitting. 
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FIG. 4. The phase diagram for d = dx, B = B^z. The shaded 
regions are topologically nontrivial which support unpaired 
Majorana bound states. In each panel the two parameters 
that are not shown take values as eq = 2 meV, = 1 Tesla, 
A|| = A-d = (to/t'o) X 1 meV, and <j> — tt/IS. Over the a range 
of a few Tesla the Zeeman energy is negligible, as seen in the 
top left panel. (We have taken the effective Bohr magneton 
to be the same as the bare value 5.7 x 10"'^ meV/Tesla.) 
The gray scale shows the value of |A| — |Imt| ~ the energy 
gap, where white is zero and black is to x 3.5 meV. The line 
indicates where the gap vanishes. 



= /iB • B, with an effective Bohr magneton /^b, and 
the spin-orbit interaction hso — \ - (cr xp) with the spin- 
orbit coupling A. A summary of the results follows with 
details given in Appendix |Bj 

The intra-layer tunneling due to the spin-orbit inter- 
action is zero since the zero energy solutions in the same 
layer have equal orbital angular momentum. The intra- 
layer Zeeman tunneling also vanishes when B _L d. This 
is true to all orders of perturbation theory, since in this 
case, the effect of the Zeeman term is to shift the chem- 
ical potentials of each spin component by ±fiBB. Since 
the two spin blocks are decoupled, we still find two zero 
modes. 

The tunneling amplitudes are 

Ret — — 2ioMBB • dcos(^ + 2iQ(d x A)^ siniyj, (13a) 
Imt = 2io£o sin (13b) 

Re A = 2tQX ■ d± cos ip + 2ioMB(B x d±)z sin ipll3c) 

Im A = 2fQdj^ • (d X A) cos tp 



+ 2ioA*B 
ifisB ■ d 



d_L X (B X d) 



sin(p. 



(13d) 
(13e) 
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Here Lp = ($i — $2)/2 is half of the superconducting 
phase difference between the layers, Eq is the energy scale 
of direct electron tunneling between the layers, and Iq 
and tg are, respectively, direct and differential overlap 
integral between the layers. We use the notation that a 
is the opposite of cr, the signs (—1)^ = 1 = —( — 1)"^, and 
dj^ is the unit vector along the projection of d onto the 
xy-planc. 

When ip — 0,we see that Ret = — 2toMBB-d, Imt = 0, 
ReA = 2t'f^X-dj_ andlniA = 2t^d_L • (d x A). There is no 
supercurrent, but Kitaev's condition Q yields B • d ^ 
and 1^0 1 > 1, which is unlikely to be satisfied since to oc 

g — a/ao 

For (f ^ 0, it becomes possible to satisfy Kitaev's con- 
dition. The supercurrent is also nonzero in this case. 
Kitaev's condition holds if 



^b\B • d| < to |(d X X)^_sm(p\ 



(14) 



which dictates that d must have nonzero components in 
the xy-plane and that A also have components normal 
to d_L in the xy-plane. The additional condition ^ for 
the supercurrent could also be satisfied if either (say, for 
small if) 

\t'f^X • d_L| > |to£ol and/or |t^d_L • (d x A)| > |to£ol, 

which means that A should have some components in the 
xy-plane, or 



/xb|B X d_L| > |eol 



(16) 



Since Zeeman energy is generally very small (^ 10^^ eV 
at 1 Tesla) the last condition would require fine tuning 
to hold. 

Our discussion suggests that unpaired Majorana 
modes may be stabilized if a supercurrent flows between 
the layers, d is unlocked from the z-axis, B • d — 0, 
and there is a spin-orbit coupling in the xy-plane. To 
illustrate let us take d = d_L to be in the xy-plane, and 
B BzZ. Then (14) is automatically satisfied if A has 



normal components to d in the xy-plane. If in addition 
A has components along d, it is possible to satisfy ( 15 ) 
for large enough spin-orbit coupling. See Fig. |4] 



IV. EXCHANGE STATISTICS OF VORTICES 



the middle with Majorana bound state 7. Using the anti- 
commutation of Majorana operators, we can see that for 
two incoming vortices U" and U' do not commute. In 
fact, they anticommute. In a vortex interference experi- 
ment, this non-A belian structure makes the interference 
term vanislP^'^B^ 

What happens in the layered case? The vortex is now 
an extended object with many Majorana bound states. 
When it encircles another vortex, the Majorana fermions 
on each layer acquire a minus sign, generated by 



N 



U'n 



(17) 



It is easy to see that and C/^ do not commute only 
when N is odd. This "odd-even" effect in the number of 
layers is a robust feature that depends on the existence 
or not of an unpaired Majorana fermion within the loop. 

Now let us consider a moving FQV with two Majorana 
bound states in each layer. First we look at the case when 
all Majorana tunneling are turned off. When the FQV 
goes around a stationary HQV with a single Majorana 
bound state per layer in, say, 7^^, it was recently showrl^ 
that j-fi — )■ —^^i while — 7^^. The generator of this 
transformation is then, f7j^ = nfci7t€7T^- Since this is 
in the same form as (17 1, the results are also the same: 



the exchange of a FQV and a HQV is non-Abelian only 
when N is odd. 

When Majorana tunneling amplitudes are turned on, 
the situation is quite different. Now, the Majorana chain 
in the FQV has two topologically distinct phases. In the 
trivial phase, all the Majorana fermions are paired into 
regular fermions that are all either occupied or unoccu- 
pied. Therefore, in this case the Majorana states become 
decoherent and non-Abelian statistics is lost over a long 
enough time. In the nontrivial phase, there are unpaired 
Majorana fermions at the ends of the vortex lines on the 
surface of the superconductor. Then, the non-Abelian 
statistics between a FQV and a HQV is recovered for 
odd N similar to the case of a single layer. An even more 
interesting, but rather speculative situation is when the 
vortex is partially pinned so that only one of the end 
Majorana fermions can move, in which case non-Abelian 
statistics is obtained at a single surface layer regardless 
of the evenness or oddness of N both between FQVs and 
between a FQV and a HQV. 



Majorana bound states in superconducting vortices act 
as non-Abelian anyons in a single layer^'^. The reason is 
that when a vortex is taken around another in a loop, 
the superconducting phase changes by 27r. Therefore the 
Majorana bound states 7', 7 at the vortices each acquire 
a minus sign. This is generated by the transformation 
U' = 7'7 on the zero-energy Hilbert space. Let us now 
consider an experiment where a beam of vortices with 
Majorana bound states {...7", 7'} (say, HQVs in a p- 
wave superconductor) circumvents a stationary vortex in 



V. EXPERIMENTAL REMARKS 

Since the unpaired Majorana states discussed in this 
work appear at the surfaces of the sample, one might ex- 
pect that surface sensitive probes such as scanning tun- 
neling spectroscopy would be able to detect them. The 
tunneling spectrum would show an enhanced resonance 
at zero bias, which exponentially decays away from the 
Majorana bound state. A possible setup is shown in 
Fig. [1] When the topological superconductor has an even 
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number of layers, there is in total always an even num- 
ber of Majorana fermions bound to the hole boundaries. 
In the trivial phase, all Majorana fermions are paired up 
and there will be no zero-bias resonance. In the nontriv- 
ial phase, however, two of these Majorana fermions are 
bound to the top and bottom layers and are separated by 
an energy gap from the rest of the spectrum. The tun- 
neling spectrum would then show a zero-bias resonance. 
As we have discussed, such a state may be obtained in a 
chiral p-wave superconductor when a supercurrent flows 
between the leads and along the hole. 

A fine energy resolution is needed to detect the zero- 
bias Majorana resonance. First, one needs to be able to 
resolve the mini-gap ^ v\/R, separating the Majorana 
state from the rest of the vortex core bound states, where 
R is the radius of the hole (for a pinned vortex, R is re- 
placed with the superconducting coherence length) . This 
mini-gap usually ranges in the sub-meV. In addition, one 
must resolve the gap that separates the unpaired Majo- 
rana bound state from the rest of the spectrum in the 
Majorana chain along the vortex line. This gap goes as 
~ |A| — |Imi| and is shown in the grayscale of Fig. |4] 
In our estimate, this gap would also range in the sub- 
meV. However, the weight of the states above this gap is 
greatly diminished since they are extended along the hole 
and therefore have a smaller overlap with the tunneling 
tip. Overall, sub-meV energy resolution and similarly 
low temperatures would be typically needed for such an 
experiment. 

These authors have recently proposed^^ an interfer- 
ence experiment based on the Aharonov-Casher effect 
to detect the non-Abelian statistics associated with un- 
paired Majorana fermions. The present work shows that 
a setup similar to that shown in Fig. [T] with unpaired 
Majorana fermions bound to the top and bottom lay- 
ers is more suitable for such an experiment. We note 
that the ability to grow high-quality thin film samples 
of Sr2Ru04 in a layer-by-layer fashion has been recently 
demonstrated^^. As our discussion shows, such a degree 
of control will be very useful in the clear identification 
of Majorana fermions and their potential applications in 
layered superconducting material. 



VI. CONCLUSION 

We studied the problem of Majorana modes on a vor- 
tex line threading a layered two-component topological 
superconductor. We mapped the problem to that of a 
general quantum wire (or a Majorana chain with two 
Majorana states per site) with nearest neighbor pairing 
term and a finite supercurrent. We established both an- 
alytically and numerically that the topological phase of 
the quantum wire, which supports unpaired Majorana 
bound states at the boundaries of the chain, survives in 
the presence of a supercurrent when the latter is below 
the critical current of the wire. This generalizes the con- 
dition derived by Kitaev for unpaired Majorana modes^ 



(that the quantum wire be partially filled). 

For a spin-triplet chiral p-wave superconductor we cal- 
culated the Majorana tunneling amplitudes due to spin- 
rotation breaking interactions and derived the conditions 
for which unpaired Majorana modes may be stabilized in 
a vortex line at the top and bottom layers. Our study 
suggests that unpaired Majorana bound states could ex- 
ist even when Majorana tunneling in and between layers 
is nonzero if a supercurrent flows between the layers and 
for certain spin-orbit couplings. These vortices satisfy 
non-Abelian braiding statistics with half-quantum vor- 
tices when there is an odd number of layers. More gen- 
erally, our results suggest that the supercurrent could be 
useful as an additional control knob for realizing unpaired 
Majorana fermions in one-dimensional systems. 
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Appendix A: Majorana fermions in a chiral p-wave 
superconductor 

1. Majorana edge modes 

For a spin-polarized Px + ipy condensate, 

A==iwAe"*/^(9:,-i9j,)e**/^ (Al) 

This condensate has a single Majorana mode in each layer 
at the boundary of the system normal to the xy-plane 
r = (r, 9). For a disk of radius i?, threaded by n vortices, 
i.e. when $ = nO^ the energy of these modes \s E = 
{vj\/R)l (see Appendix [A| for details) where the angular 
momentum I E Z for odd n and I G Z + ^ for even n. 
Thus, for a single vortex there is a single zero energy edge 
state for any R. 
We take 

h = -^V2 + -^pI + K(r) + yp(z) - A^e, (A2) 

with the confining potential (r) chosen for a disk of ra- 
dius R to be Vc{r <^ R) and Vc(r > i?) = where 
Vo is large and positive. The periodic potential Vp{z) cre- 
ates the layered structure: the layers are potential min- 
ima ai z = ia where i is an integer. We first concentrate 
on a single layer, i.e. we solve \p1/2m + Vp{z)\ ^/Jzj{z) — 
Ez^z,e{z) around each z = la and project to the lowest 
energy state. Next, we define the electrochemical poten- 
tial e(r) = fic — Voir) — E^, and choose e{R) = so that 
e will be negative on one side of the edge (the insula- 
tor at r > R) and positive on the other (the topological 
superconductor at r < R). Let us specialize to the spin- 
polarized chiral p-wave condensate in (|A1|. Linearizing 
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the Hamiltonian near r = R, taking the phase of the or- 
der parameter to be $ = n9 where n is the number of vor- 
tices through the superconductor, and choosing the an- 
gular dependence of the spinor to be exp [i(Jz{n — 1)6/2], 
we can factorize the BdG Hamiltonian as 



e(r) -I- v^Ty dr + 



1 

2r 



H 



The radial part contains 
[e(r) - VA {dr + ^)] /o = 0, 



, 1 

Jo = ^ CXp 



1 



e{p)dp 



(A3) 
mode, 

(A4) 



We assume that the potential is sharp enough so that 
higher energy states are not important. We can now 
write the Nambu spinor projected into this zero energy 
state 



i)^/2e-/Vc(r)^(0) 



1 



(A5) 



where = $ — 716' is a constant. The effective Hamilto- 
nian for the Hermitian field i!{9) is = —i{vA/ R)dg. 
We can expand ip{0) = e^'^V'ij where ipi satisfies tpj — 
ip-i. The energy levels are Ei = —{va/R)1- We note that 
the spinors are single- valued only if Z ± (n — l)/2 € Z. 
Consequently, if n is odd, then / € Z, and if n is even, 
l£Z + 1/2. 



As is apparent from equation (A3), the eigenvalue of 



(Ty is chosen by the asymptotic behavior of the electro- 
chemical potential e, locking the direction of propagation 
of the edge to the spatial profile of sgn(e(r)). 



2. Majorana bound states 

For the HQV in the spin-triplet chiral p-wave conden- 
sate in (12), we take d{r,9) — (cos |, — sin |, O) , and 
<f> = + 9/2 with a constant 4>. Since 



e*^/2(cr.d)ifT,y 



(A6) 



only the spin-up component sees a vortex winding in the 
phase. The BdG equation has a zero mode 



fir) = Ue 



e < ^mv 



A; 
2 

A' 



(A7) 
(A8) 



where ~ (1,0)-^, k = ^^2m\e — |mw^| and A/" is a 
normalization factor. Here, we also assumed a vanishing 
small vortex core so that va is constant everywhere but 
at the origin, as well as a constant e away from the edges. 



The d vector in the presence of a FQV does not con- 
tain any winding and we simply take it to be a constant. 
The Hamiltonian is invariant under a SU(2) spin rota- 
s' 



tion * U^', with U 



S* 



S e SU(2), and 



d Rg d', where Rs € S0(3) is the rotation corre- 
sponding to S. With this symmetry we can rotate d 1— > y, 
for which the Hamiltonian decouples into spin compo- 
nents with a vortex phase winding ^ — (j) + in both 
spin components. There are two zero modes 
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/ 



-iVn 



(A9) 



with (7 = t,i and vi = (0,1)"^. Rotating back to the 
original basis, the full solutions are found as Xa = Ux'a- 

Appendix B: Majorana tunneling 

1. Intra-layer Majorana tunneling 

The BdG Hamiltonian receives a contribution iJz from 
Zeeman interaction, which after the rotation Rs : d i— > y 
reads 



U^HzU 



fJ'B 







-/iB(5ta.B5r ' • ^^^^ 



So, using the identity 

5^ • BS = o- • i?sB 

we have 



a = -z(x^|C/ti/zC/|4) 



-2/xbB • d. 



(B2) 



(B3) 



The spin-orbit term when A||z was considered in Ref.!^ 
and found not to lift the degeneracy of the two Majorana 
modes. We show that, to the first order in perturbation 
theory, this is actually true for general A. We have 



_ , A - (cr X p) 
- ' -A • X p)^ 



(B4) 



Therefore, the overlap (xtl^solx^) "^iH contain terms 
involving components of the integral J{fpf)dr (x 
/(V|/p)dr = since / vanishes at infinity. We conclude 



{xt\Hso\xi)=0. 



(B5) 



When B • d — 0, the Zeeman tunneling vanishes, 
which is true to all orders of perturbation theory. The 
chemical potential is shifted for each spin component to 
fij. — fi^ ± fi^B. Since the two spin blocks are decou- 
pled, we still find two zero modes but with new functions 
f±. If we then turn on the spin-orbit interaction, to the 
first order of perturbation theory, / /+(A x 'V)xf-dr = 
since / f+dzf-dr — and J f^dyf^dr = by cylindrical 
symmetry of f±. So (|B5| still holds. 
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we find (to) erg- = and 



(to)crcr = toSoSillip, 



(B7) 



where tp — ($i — $2)/2 is half of the superconducting 
phase difference between the layers. So, in a uniform 
condensate to = 0, even though Cooper pairs move freely 
in all three direction! 

The Zeeman interaction yields the tunneling ampli- 
tudes 



FIG. 5. (color online) The geometry of the d vector. 



2. Interlayer Majorana tunneling 

Let us consider the overlap between the Majorana 
modes of a FQV in two adjacent layers indexed as 1 
and 2. It is given by t = —i {x2\Ho\xi) i where Hq = 
I 



toAifiB -Im {e'^^v^S^crSv^,) 



(B8) 



-/ 



with eg energy scale for the tunneling 



process. So, by defining 



to = 2y" /2(r)/i(r)drdz, 



(B6) 



We shall use the following parameterization d = 
(cos (5 sin /3, cos (5 cos /3, sin 5) visuahzed in Fig. [s] and de- 
fine d_L — (sin/3, cos/3, 0) as the unit vector along 
the projection of d onto the xy-plane. Then S = 
exp{i/3az/'2) exjp{—i5ax/'2). It is straightforward to find 



{tz)aa = {—^ytofiB{—Bx smSsiu (3 — By sin (5 cos P + cos 6) sinip (B9) 
{tz)aa = —{—^ytof^B{Bx COS (5 sin /3 -I- By cos (5 cos /3 -I- B^ sin<5) cos (p + tofisiBx cos/3 — By sin /3) sin ip, (BIO) 



Substituting back the components of d and dj^ we have then 

{tz),, = (-l)"toMB \di_ X (B X d) 



sm ip, 



(Bll) 
(B12) 



(tz)aw = *oMb (B X d_L)z sin ip - (-l)'^B ■ dcos ip 
The contribution from the spin-orbit interaction can be similarly written as 

{tso),,,=-l{X2a\U^HMla') 

= vl [e'fS^trS + e-''^{S'^(TSf] v„. ■ J /2(A x V)/idrdz. 

= 2Re {e^^v'^S^'crSv^,) ■ J /2(A x V)hdrdz. 

Assuming cylindrical symmetry (i.e. that the vortex does not wiggle), we have / f2{dx,dy)fidr — 0. Then (B13) 
simplifies to 

{tso)aa = (-l)'^to(Aso,xSin(5cos/3- Aso,i;Sin(5sin/3)sini^ (B14) 
sin^ -I- Aso,j, cos/3) cos p + (— l)'^to(— Aso,^ cos 5 cos (3 + Xso,y cos (5 sin /3) sinp, (B15) 



(B13) 



with 



Then 



t'o^^J f2djidrdz. 



(tso),, = (-l)%di-(dxA)cosv?, 
{tso)^^^to A-dj_cos.^ + (-l)°'(dx A)^sin.^ 



(B16) 

(B17) 
(BIB) 
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